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H Abstract 



Recently, a color decomposition using structure constants was introduced for purely gluonic 
tree amplitudes, in a compact form involving only the linearly independent subamplitudes. We 
give two proofs that this decomposition holds for an arbitrary number of gluons. We also present 
and prove similar decompositions at one loop, both for pure gluon amplitudes and for amplitudes 
with an external quark-antiquark pair. 



Submitted to Nuclear Physics B 



1 Address after November 1, 1999: Theory Division, CERN, CH 1211 Geneve 23, Switzerland 
2 Research supported by the US Department of Energy under grant DE-AC03-76SF00515 



1 Introduction 

The computation of multi-parton scattering amplitudes in QCD is essential for quantitative pre- 
dictions of multi-jet cross-sections at high-energy colliders. The Feynman diagrams for such an 
amplitude can generate a thicket of complicated color algebra, tangled together with expressions 
composed of kinematic invariants. An extremely useful way to disentangle the color and kinematic 
factors is via color decompositions ||, [21 1| Q, || in terms of Chan-Paton factors M, or traces of SU(N C ) 
matrices in the fundamental representation, A a . The traces are multiplied by purely kinematical co- 
efficients, called partial amplitudes or subamplitudes. Originally motivated by the representation of 
gluon amplitudes in open string theory, such decompositions have been widely studied and applied 
at both the tree level ||| and the loop level 0. Two major advantages of the trace-based color de- 
compositions are that (a) subamplitudes are gauge-invariant, and (b) an important 'color-ordered' 
class of them receive contributions only from diagrams with a particular cyclic ordering of the ex- 
ternal partons; these subamplitudes therefore have much simpler kinematic properties than the full 
amplitude. 

Although the standard color decompositions are very effective, in some cases they are not quite 
optimal. For example, the decomposition of the purely gluonic tree amplitude is overcomplete. 
Explicitly, the n-gluon amplitude is written in terms of [n — 1)! single-trace color structures, as 

AT(9i,92,...,gn)=9 n - 2 E tr(X^X^-..X a ^)AT(l,a 2 ,...,a n ), (1.1) 

fTgSn-1 

where A^ ee are the subamplitudes. The cyclic invariance of the trace has been used to fix its first 
entry, and the sum is over the set of non-cyclic permutations of n elements, a £ S n -\ = 5 n /Z n , 
where <r(2, . . . , n) = (02 > • • • , cr n ). However, the (n — 1)! subamplitudes appearing in the equation are 
not all independent. Besides cyclic and reflection invariances, which they inherit from the traces, 
subamplitudes also obey a U(l) decoupling identity (also called a dual Ward identity, or subcyclic 
identity) |, |, 

< ce (l,2,3,...,n)+< ee (l,3,2,...,n) + --- + < ee (l,3,...,n,2) = 0. (1.2) 

This identity can be derived by setting A° 2 equal to the unit matrix and collecting terms containing 
the same trace; this gives the amplitude for a U(l) 'photon' and n — 1 SU(N C ) gluons, which must 
vanish. More general identities can be derived by assigning the SU(N C ) generators for the external 
gluons to commuting subgroups such as SU(Ni) x SU(N2) ||. Kleiss and Kuijf found a linear 
relation between subamplitudes |1(| (see section ||), which is consistent with all of these identities, 
and which reduces the number of linearly independent subamplitudes from (n — 1)! to (n — 2)!. 

In contrast, the conventional color factors for gluonic Feynman diagrams are composed of struc- 
ture constants f abc . In this case, the U(l) and generalized decoupling identities are all manifest; 
for example, the U(l) generator commutes with all other generators, so all structure constants con- 
taining it are zero. On the other hand, a particular string of contracted structure constants will 
not appear ab initio with a gauge-invariant kinematic coefficient. To pass from the / -based de- 
composition to the trace-based decomposition, one substitutes f abc = — itr (A a A b A c — A fc A a A c ) and 
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simplifies and collects the various traces. In the process, gauge invariance of the partial amplitudes 
is gained, but properties such as U(l) decoupling are no longer manifest. Another disadvantage of 
trace-based decompositions at tree level is that the color structure of amplitudes in the multi-Regge 
kinematics is obscured f|Tl|. 



The trace-based decomposition is also not optimal at the loop level. At one loop, the standard 
color decomposition || for n-gluon amplitudes includes double trace structures, tr(A Ql • • • A" 0-1 ) 
x tr(A° c • • • A a "), in addition to single traces of the type that appear at tree level. The subamplitudes 
multiplying the single trace structures give the leading contributions in the large N c limit, and they 
are color-ordered. The double-trace subamplitudes, corresponding to subleading-in-iV c contributions, 



can be written in terms of permutations of the leading (color-ordered) subamplitudes [12] in a formula 
reminiscent of the tree-level Kleiss-Kuijf relation (see section |||). Because the permutation formula is 
rather complicated, its implementation in a numerical program can be rather slow. Similar formulae 
hold for one-loop amplitudes with two external quarks and n — 2 gluons [jl3[ . 

Recently, a new color decomposition for the n-gluon tree amplitude, in terms of structure con- 
stants rather than traces, has been presented [j~4j| , whose kinematic coefficients are just color-ordered 
subamplitudes. Because of this property, the new decomposition retains all the advantages of trace- 
based decompositions, yet avoids the disadvantages mentioned above. In particular, the tree-level 
n-gluon decomposition is automatically given in terms of the (n — 2)! independent subamplitudes, 
in a form that is very convenient for analyzing the high-energy limit. 

The purpose of this paper is to derive the new tree-level decomposition in two different ways, 
and to present similar f a -based color decompositions at one loop, whose kinematic coefficients are 
again color-ordered subamplitudes. The leading and subleading-in-iV c contributions combine neatly 
into one expression in the new decompositions. Also, gluons circulating in the loop are put on the 
same footing as fermions in the loop. 

Parton-level cross-sections require amplitude squares or interferences which are summed over all 
external color labels. We provide general expressions for color-summed cross-sections in terms of the 
new color decompositions. In appendix [A] we give explicit evaluations for quantities encountered in 
cross-section computations through 0{cQ). 

The paper is organized as follows. In section |2| we give two proofs of the new color decomposition 
for the n-gluon tree amplitude conjectured in ref. [O], and provide accordingly the square of the 
tree amplitude summed over colors. In section || we present (and prove) new color decompositions 
for one-loop n-gluon amplitudes and one-loop amplitudes with an external quark-antiquark pair 
plus n — 2 gluons. In addition, we compute the color-summed interference terms between tree 
amplitudes and one-loop amplitudes, and the square of one-loop amplitudes, which are relevant 
respectively for next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) calculations 
of jet production rates. The color-summed interference terms and squares are given in terms of 
color matrices; appendix A provides many of those required for jet rate computations up to 0(a\ ). 
Finally, we outline how a new color decomposition for one-loop amplitudes with external photons, 
gluons and a quark-antiquark pair can be obtained from the amplitudes with only gluons and a qq 
pair. In section |] we summarize our present understanding of the color decomposition of tree-level 



and one- loop amplitudes, and comment briefly on possible extensions to multi-loop amplitudes. 

2 Tree Color Decomposition 

The new color decomposition for the n-gluon tree amplitude is[| [14] 

AT e (gi,---,g n ) = (ig) n - 2 £ f**>^f**°*** ■■■ r-^-i-AT e (i,a 2 ,...,a n . l ,n), 

erG<S n -2 

= g n - 2 J2 (F a ^---F a ^)a 1 a n A t ^(l,a 2 ,...,a n „ 1 ,n), (2.1) 

(TGSn-2 

where {F a )i, c = if bac is an SU(N C ) generator in the adjoint representation. This color decomposition 
is analogous to the standard decomposition for the tree amplitude with a quark-antiquark pair and 
n-2gluons @|, 

A t 7 r(q 1 ,g 2 ,...,g n - 1 ,q n ) = g n - 2 £ (\°">* ■ ■ • A a ^-i)>A^ e (l g ,a 2 ,... ,<7 n _i,n f ) . (2.2) 

(TG5 rl _2 

The only difference between the two is the representation used for the color matrices, namely the ad- 
joint representation for the n-gluon amplitude and the fundamental representation for the amplitude 
containing a qq pair. 



Eq. (2.1) is written in terms of the (n — 2)! subamplitudes where legs 1 and n are adjacent. This 
is precisely the basis of linearly independent subamplitudes which is singled out by the Kleiss-Kuijf 
relation. In fact, we shall show that eq. ( [2.1| ) is equivalent to the Kleiss-Kuijf relation JO]], which 
can be written as, 

4T(l,{a},n,{/3}) = (-1)»* £ 4T(l,a({a}{/? T }),n). (2.3) 

aeOP{a}{/3T} 

Here {a} U {(3} = {2,3, . . . ,n — 1}, np is the number of elements in the set {/?}, the set {/3 T } is 
{/3} with the ordering reversed, and OP{a}{/3 } is the set of 'ordered permutations' (also called 
mergings |1(J) °f the n — 2 elements of {a} U {(3 } that preserve the ordering of the a« within {a} 
and of the /3j within {/3 T }, while allowing for all possible relative orderings of the «j with respect to 
the Pi. 

We wish to show that the new color decomposition (^J) is equivalent to inserting the Kleiss- 



Kuijf relation (2.3) into the standard color decomposition (|1.1| ). We first substitute for the structure 
constants appearing in eq. (2.1), f abc = — itr(A a [A 6 , A c ]), and use the identity 



fttia 2 ii M10312 . . . fin-3«n-i«n = (— t) n— 2 tr (A ai \\ a2 \\ az [A a ™ _1 A a "l 11) (2 4) 

We want to identify all terms that contain a trace of the form 

tr(A ai A QQ i ...\ a ^-2- q \^ x a Pl ...x a 0i), (2.5) 

1 We choose the normalization of the fundamental representation matrices as tr(A a A i> ) = S ab . Hence our structure 
constants f abc are larger than the conventional ones by a factor of y/2. 



because these should give rise to j4^ ree (l, {a}, n, {/3}). Since there are n — 2 commutators, there are 



2™ terms on the right-hand side of eq. ( |2.4D , but only ( ~ ) of them have exactly q A matrices 
appearing to the right of A a ™. The aj indices on the matrices to the right of X a " must come 
in reversed order compared to how they appear in the f abc string, but they can appear in that 
string in any relative order with respect to the a^ that end up to the left of A a ™. Thus, for any 
ordered permutation a G OP{a}{/3 T }, the subamplitude A t ^ ee (l,a({a}{/3 T }),n) appears in the 
new decomposition ( |2,lD accompanied by the desired trace (|2.5| ); a relative sign of (— l) q = (—l) n ' s 
comes from the commutators. Collecting all such ordered permutations, we obtain the Kleiss- 
Kuijf formula for ^ ree (l, {a}, n, {/?}), eq. (\L3j), thus establishing its equivalence to the new color 



decomposition (l2.l|) . 

Hence we can derive the new n-gluon color decomposition via its connection to the Kleiss-Kuijf 
relation. The latter relation was checked up to n = 8 in ref. [flOf] . It can be proved for all n using 
the same techniques that were used to prove an analogous one-loop formula fT2| , eq. ( JJT2J ) below. 
Consider first the set of color-ordered Feynman diagrams with only three-point vertices, which are 
'multi-peripheral' with respect to 1 and n, by which we mean that all other external legs connect 
directly to the line extending from 1 to n; i.e. there are no non-trivial trees branching off of this 
line. Label the legs on one side of the 1-n line by on, and those on the other side by (3j. These 
diagrams contribute to both sides of the Kleiss-Kuijf relation ( |2.3| ) in the proper way: There is no 
relative ordering requirement on the on with respect to the j3j on the left-hand side of the relation, 
and the ordering is summed over on the right-hand side; the sign factor (— l) n ' 3 comes from the 
antisymmetry of color-ordered three-point vertices (in Feynman gauge). 

Next consider diagrams with non-trivial trees attached to the 1-n line. They work in exactly 
the same way as the multi-peripheral diagrams if the leaves (external legs) of each tree all belong 
to the same set, either {a} or {/?}. However, if a tree contains leaves from both sets, then the 
diagram does not appear on the left-hand side of eq. ( |2.3|) , so one must show that it cancels out 
of the permutation sum on the right-hand side. For trees with only three-point vertices, this can 
be done using the antisymmetry of the vertices. For other cases the cancellations are slightly more 
complicated; one way of establishing them is via the generalized decoupling identities ||] mentioned 
in the introduction III 
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Figure 1: Graphical representation of a multi-peripheral color factor. A vertex stands for f abc , and a bond 
for S ab . 



There is another way to prove eq. fl2.1|) . For this argument it is convenient to use a graphical 
notation for color factors made out of structure constants, in which f abc is represented by a three- 
vertex and 5 (an index contraction) is represented by a line ||]. Then the color factors appearing 



in eq. ( |2.1[ ) are associated with the multi-peripheral color diagrams shown in fig. |l[ The color factor 
for a generic n-gluon Feynman diagram is not of this form. However, we can use the Jacobi identity 



rdac rcbe 



rdbc fcae 



rate rdce 



(2.6) 



shown graphically in fig. g, to put it into this form ]l]]. Fig. y represents the color factor for a generic 
Feynman diagram graphically, as a tree structure with only three-point vertices. It also shows how 
it can be simplified into multi-peripheral form by repeated application of the Jacobi identity. (If 
the line running from 1 to n is in the fundamental representation, then the same steps, but using 
A a A 6 — A fe A a = if abc \ c , lead directly to the color decomposition $L% ) for a quark-antiquark pair and 
n — 2 gluons.) 



Figure 2: The Jacobi identity for structure constants in graphical notation. 




n 




Figure 3: Graphical representation of the color factor for a generic tree-level Feynman diagram. Also shown 
is a step in its conversion to multi-peripheral form, by using the Jacobi identity in the dashed region. 



So far we have established that a color decomposition of the form (f2T|) exists, but we have not 
yet shown that the kinematic coefficients are equal to the partial amplitudes A^ ee (l, 02, . . . , cr n _i, n). 
However, we can do this rather simply by equating the new decomposition (but with unknown kine- 
matic coefficients) to the standard one, eq. ( |l.ip , and contracting both sides with tr(A ai A a<T 2 • • • \ a(y n-i A an ), 
i.e. summing over all cij, i = 1, ... ,71, for some permutation a £ S n -2- Because the kinematic co- 
efficients contain no N c dependence, we may retain only the leading contraction terms as N c — ► 00. 
But it is well known that single traces for different cyclic orderings are orthogonal at large N c Q. 
Similarly, fig. |I| shows that only one of the multi-peripheral color factors in eq. (2.1) survives the 
contraction, because the other S n _2 permutations give rise to nonplanar, and hence A r c -suppressed, 
color topologies. Thus the contraction selects a unique term from either side of the equation, cor- 
responding to the same ordering of the color indices dj, and with the same weight at large N c , 



namely JV". This proves that the coefficients in the new color decomposition (2.1) are indeed the 
subamplitudes A^ ree (l, 02, • • • , cr n -i,n). 
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Figure 4: Contraction of different multi-peripheral color factors with a single A trace. The heavy line 
with an arrow denotes the fundamental representation. The left diagram shows the contraction where aj, 
i = 2, 3, . . . ,n — 1, appear in the same order in the A trace as in the multi-peripheral color factor. Every other 
contraction gives rise to a nonplanar (and hence color-suppressed) diagram of the type shown on the right. 



Because we showed that the new color decomposition is equivalent to the Kleiss-Kuijf relation, 
our second proof of the new decomposition can alternatively be viewed as a proof of the Kleiss-Kuijf 
relation. Interestingly, this proof does not rely on any kinematic properties of color-ordered Feynman 
diagrams; only group-theoretic properties such as the Jacobi identity and the large-iV c behavior of 
color traces were used. 

The decomposition (^Tj) has a color-ladder structure which naturally arises in the configurations 
where the gluons are strongly ordered in rapidity, i.e. in the multi-Regge kinematics ]15[. (Indeed, 
this is how eq. (2.1) was first discovered IJUfl .) In contrast, some laborious work is necessary to 
obtain the color-ladder structure from the trace-based decomposition ( |1.1| ) in the multi-Regge kine- 
matics [11 1 . Let gluons 2, . . . , n— 1 be the particles produced in a scattering where gluons 1 and n are 
the incoming particles. The available final-state phase space may be divided into (re — 2)! simplices, 
in each of which the produced gluons are ordered in rapidity. Let us take the simplex defined by the 
rapidity ordering yi > 1/3 > • • • > y n ~2 > 2/n-i, and consider its sub-simplex with strong rapidity 
ordering j/2 3> V3 S> ■■■ S> J/n-2 S> Vn-i- Then the crucial point to note is that the dominant 
subamplitudes in eq. ( |L1| ) in the multi-Regge limit are A^ ee (l, {a}, re, {/? T }), where {a} and {/3} 
are both increasing sequences, whose union is {2,3, . . . ,re — 1} ({/3 } is {/?} in reversed order). In 
other words, OP{a}{/3} contains the sequence {2, 3, . . . , re — 1}. If q is the number of elements in 
{/3} , then for each q there are ( n ~ ) such choices of {a}, {/3}. Summing over q = 0, 1, 2, . . . , n — 2, 
there are 2 n ~ 2 orderings in total Jig] • Using the identity (|2.4| ) and the fact that in the multi-Regge 
kinematics 

< rce (1, {a}, re, {/3 T }) = (-l)Mf (1, 2, . . . , n), (2.7) 

one obtains the color-ladder structure of eq. Q2.1| ), with just one allowed string of structure constants 
/ . The procedure can be repeated for all of the (re — 2)! simplices, thus generating the (re — 2)! 
strings of f abc, s of eq. (2.1). 

The tree- level partonic scattering cross-section is given by the square of the tree amplitude A^ ee , 
summed over all colors. This expression can be written either of two ways, 



£KT(l,...,re) 



(9 



colors 






r-A 



tree 



(Af 



(2- 



(9 



(n-2)l 
2\n— 2 \ "* ~ /(tree/ /itree\* 
/ 4 c ij A i \ A j ) > 



(2.9) 



where eq. ( |2.8D is obtained from eq. (1.1), while eq. (2.9) is obtained from eq. (2.1), and the subscript 
i on vl* ree now refers to the subamplitude A^ ee evaluated for the i th permutation P, in S n —i or S n —2, 
respectively. In eq. fl2.8|) , the color matrix Cij is 

(2.10) 



Cij = ^ tr(X ai P l {X a2 ■ ■ ■ A a ™}) [tr(\ ai P j {\ a2 ■ ■ ■ A a ™})]* , 



colors 



whereas in eq. ( |2.9j ) the matrix is 

c l3 = ^ (P t {F a > ■ ■ ■ F^}) aian U Pj {F» ■ ■ ■ F^}) aian 



(2.11) 



colors 



In refs. |1C, 17], the Kleiss-Kuijf relation was used to calculate dij from c^. From eq. ( p. 11 ) we 
see that it can be calculated directly in terms of structure constants. In appendix [A] we give c^ for 
n = 4,5,6. 

Although the expression (2.9) seems to have fewer terms than eq. fl2.g|), the sparseness of the 
matrices is also important for determining which expression has the fastest numerical evaluation 17]. 
For an evaluation that is good only to leading order in N c — the Leading Color Approximation, or 
LCA — it is best to use eq. ( [2.8D and the large-iV c orthogonality of different single traces, to obtain 

1 



53i^r(i,...,n)p 



tri- 



colors 



T- 2 e n (N c ) 53 

ff€S„_i 



A tree (l,a 2 ,...,a n )\ 2 + 



7V-2 



where 



C n {N c )=N?-\N 2 c -l). 



2.12 



(2.13) 



Up to n = 5, the LCA is exact, and so eq. fl2.8| ) is also superior to eq. (2.9) for an exact evaluation. For 
n > 6, where the LCA has subleading corrections, it is best to take the leading terms from eq. ( |2.8| ), 
and the subleading terms from eq. ( |2.9| ) (these are the terms proportional to a in eq. ( |A.3| )) Jl7[ . For 
n = 7, ref. [|l^] employed the linear dependences in the Kleiss-Kuijf relation to find an even more 
compact form for the subleading terms than using eq. (2.9). (See ref. [18] for another approach to 
numerically evaluating multi-parton amplitudes beyond the LCA.) 



3 One-Loop Color Decompositions 

The standard color decomposition for one-loop n-gluon amplitudes in SU(N C ) gauge theory with nf 
flavors of quarks is |5| 



4l-loop _ r 



N c 53 tv(X a ^---X a ^)A l ^ l (a 1 ,...,a n ) 

crGSn/^n 

|n/2j+l 

+ H E tr(A^---A^-i)tr(A a ----A a -)A n;c (cT 1 ,...,(T n ) 

C=2 cr£:Sn/Sn;c 
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+ nf Y, tr(A^---A a -)4 1 . / 1 2l ( CTl ,...,a, 

cr£S„/7i„ 



(3.1) 



where A n - C are the subamplitudes, S n]C is the subset of S n that leaves the double trace structure 
invariant, and [x\ is the greatest integer less than or equal to x. The superscript [j] denotes the 
spin of the particle circulating in the loop. The subamplitudes Ajz;-, are color-ordered, and give rise 
to the leading contributions to the cross-section in the large- iV c limit. 

The subleading subamplitudes A n - C> i can be obtained from the leading ones A nl through the 
permutation sum |]l 

A n;c>l (l,2,...,c-l;c,c + l,...,n) = (-l)^ 1 £ A^a^. . . ,a n ) , (3.2) 

CT ecop{o}{/3} 

where a« G {a} = {c — 1, c — 2, . . . , 2, 1}, /3j 6 {/3} = {c, c + 1, . . . , n — 1, n}, and COP{a}{/?} is the 
set of all permutations of {1,2,... , n} with n held fixed that preserve the cyclic ordering of the ckj 
within {a} and of the /3j within {/?}, while allowing for all possible relative orderings of the c^ with 
respect to the /3j. 

In a dual representation of one-loop amplitudes, the subleading subamplitudes A n - c> i({a T }; {/3}) 
come from the annulus diagram where gluons belonging to {a} and {/?} are emitted from the inner 
and outer boundaries of the annulus, respectively. Thus their color properties are very similar to 
those of the tree subamplitudes A^ ee (l, {a},n, {/3}) if legs 1 and n are 'sewn together', i.e. contracted 
with S aian . Indeed, with this interpretation, eq. ( |3.2| ) has a very similar structure to the Kleiss-Kuijf 
relation ( |2.3| ). These remarks suggest that the new tree- level color decomposition should have a 
one-loop analog, where the strings of structure constants that appear are 'ring diagrams' instead of 
multi-peripheral diagrams, as depicted in fig. ||a. 




(a) -4 V- (b) 




Figure 5: (a) Ring color factors for one-loop n-gluon amplitudes, (b) The corresponding color factors for 
one-loop amplitudes with an external quark-antiquark pair. 

The new one-loop n-gluon color decomposition, expressed in terms of adjoint generator matrices 
(F a ) bc , is 

A' lo ° P = 9 n E [ tr ( Fa<T1 • • • paa " ) A Si fa, • • • , *•„) + 2n/ tr(A a -i • • • A°*» )A% 2] (a 1 ,...,a n )] (3.3) 



Here S'n-i = S n /Z n is the group of non-cyclic permutations, and 72 is the reflection: 72.(1, 2, . . . , n) = 
(n, ... ,2,1). Thus the number of linearly independent subamplitudes is (n — l)!/2. 

The terms proportional to nj in eq. (3.3) follow from eq. (|3.1|) using only the reflection identity 



satisfied by the A n -i: A^'.^n, . . . , 2, 1) = (— l) n A^.^l, 2, . . . , n). To prove that the remaining terms 



in eq. fl3.3j ) are correct, we follow closely the second proof of eq. (|2.l|) . Using the Jacobi identity, 
the generic one-loop / color factor can be turned into ring diagrams. The manipulations are 
exactly the same as those shown in fig. |3[ once legs 1 and n are joined together to form a loop. This 
establishes eq. ( |3.3D , but with an unknown coefficient for tr(i ?<v i ■■■F a ' yri ). Next we equate this 
expression to the standard color decomposition ( |3.1| ), and contract both sides with tr(A a<T i • • • A a<Tn ). 
Again a unique term survives on each side of the equation, at leading order in N c , and this establishes 
that the coefficient is correctly given by Aj^ai, . . . , a n ). (If we had not removed the reflections 72 



from eq. (3.3), then both a ring diagram and its reflection would have contributed to the contraction; 
this would have led to a factor of \ in the coefficient of tr(F a<T i • • • F a<Tn ).) 

Eq. ( |3.3j ) puts the gluon and fermion loops manifestly on the same footing. The contributions 
of the subleading-color subamplitudes A n]C> i are neatly packaged into the ring-diagram color struc- 
tures. In addition, all of the generalized U(l) decoupling identities are automatically incorporated. 
Finally, eq. ( |3.3[ ) shows explicitly the difference between the tree-level and the one-loop color decom- 
position of the gluon sector for any number of external legs. For instance, if we take the amplitude 
for a Higgs boson plus three gluons, both eq. ( |l.l|) and the gluon-loop contribution to eq. (|3.l| ) have 
as a color factor tr(A a<T i A a<T 2 A a<T 3); the actual difference in the color structure is concealed in the 
different properties of the subamplitudes A^ ee and A^: v In contrast, eq. (|2.l| ) and eq. ( |3.3j ) show 
explicitly the difference in color structure. 

One-loop amplitudes contribute to NLO QCD cross-sections through their interference with 
tree amplitudes. Carrying out the color-summed NLO interference terms for amplitudes color- 



decomposed as in eqs. (2.1) and (|3.3|) is straightforward: 



E A(l,...,n)[A(l,...,n)]*\ NLO 



colors 



^l-loop)*\ ( 3 _ 4 ) 



2 J^ Re (A tTC ' 

colors 

(n-2)! (n-l)!/2 
2(g 2 r- 1 Re £ £ A**°\c ij {Afy + 2n f d ij {Af 2 \ 



where 



%' = E ^ Fa2 ■ ■ ■ i?a "" 1 })a 1 a„ [* (F*Pj{F°* ■ ■ ■ F^})]* , 

colors 

dij = E^i Fa2 --- i?a "" 1 »a 1 a„NA ai P,{A a2 ---A a "}r, (3.5) 



colors 



with Pi the i permutation in S n -2 and Pj the j* permutation in S n ~\/lZ. 
NNLO production rates include the virtual contributions, 

E ^(l,...,n)[^(l,...,n)]*| NNLO = E [2Re(^ trce (^ 2 - loop )*)+|^ 1 - loop | 2 ] . (3.6) 



colors colors 



For the square of one- loop n-gluon amplitudes on the right-hand side of eq. ( |3.6| ) , we obtain 

(n-l)!/2 r 



^^l-loopp = {g 2 )n £ 
colors *J=1 



CijA^ (A [ l 1 )* + 4 n f Re d^A^ (A 



+ An]c t] 4' 2 \A^Y 



(3.7) 



where 



^ = ^ tr (F ai Pi{F a2 ■ ■ ■ F a "}) [tr (F ai Pj{F a2 ■ ■ ■ F a "})]* 

colors 

dy = ^ tr (F ai P{F a2 • • • F a "}) [tr(A ai P,{A a2 • • • A a "})]* , 



(3.8) 



colors 



Cij is given in eq. ( |2.10| ), and Pj is the i th permutation in S n -\/lZ. In appendix [A| we give the explicit 
values for the color matrices Cy , dij , c^ , dij and cy that are required for cross-section computations 
up to 0(a,g). 

A similar analysis can be applied to one-loop amplitudes with an external quark-antiquark pair 
plus n — 2 gluons. These amplitudes have the standard color decomposition [13j, 



n-l 



A 1 n - loop (qi,q 2 ,g3,...,gn)=g n Y J E Gr !3W-- ,a n )A n . j (lq,2 q ;a 3 ,...,a n ), (3.9) 

j = l creS n -2/S n -j 

where the color structures Gr^ are defined by 

1m) 



Gr^(3,...,n) = N c (\ a * ■ ■ ■ A "),* 1 
Grg(3;4,...,n) = 0, 



M) 



Gr^(3,...,j + l;j + 2,...,n) = tr(A a3 • • • A a ^ ) (A^+ 2 • • • A a ")^ , j = 3, . . . ,n - 2, 



Grg* ) _ 1 (3,...,n) = tr(A a * • • • A a "H ; 



(3.10) 



and S 



n;j — ^j-i 



is the subgroup of S n _2 that leaves Gr^ invariant. 



The leading-color subamplitudes in eq. (3.9) can be expressed in terms of color-ordered ('primi- 
tive') amplitudes as 



A n - } l{i.q, 2q,6, . . . ,n) — A n ' (lq, 2q, 6, . . . , n) 2 /L n ' [lq, 2q,6, . . . ,Tl) 



' at n V 1 "?' A q> °i ■ ■ ■ i 



AT, 



n 



(3-11) 



whereas the subleading contributions are given by the permutation formula [13|, 
A n -j{lq,2 q ;3, ... ,j + l;j + 2J + 3, ...,ri) 



(-IF" 1 E A n' [1] (*(1* 2„ 3, . . . , »)) - ^A^W (a(1 _, 2 „ 3, . . . , n)) 

aeCOP{a}{f3} L 



(3.12) 
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where {a} = {j+1, j, . . . ,4,3}, {/?} = {1, 2, j + 2, j+3, . . . ,n-l, ra}, with 1 ? held fixed. Ineqs. ( pl|) 
and Q3.12Q , the quantities A n and A n are defined in terms of color-ordered Feynman diagrams, 
much like the n-gluon subamplitudes A„.i. However, because an external qq pair is present now, one 
also has to specify which way the line running from the antiquark to the quark turns when it enters 
the loop, left or right; this accounts for the additional index, L or R. Again the index [j] denotes the 
spin j of a particle circulating around the loop, for those graphs in the primitive amplitude where 
the external qq line does not enter the loop. The reflection identity, 

A^\l q ,3,A,...,2 q ,...,n-l,n) = (-l) n A^(l q ,n,n-l,...,2 q ,...,A,3), (3.13) 

relates the left and right subamplitudes. 

There are two types of subleading-color contributions in eq. ( 3.12J ), distinguished by whether they 



contain a factor of nt or not. Although both types are generated by cyclicly ordered permutations, 
their origins are quite different. In particular, the subleading nj terms arise from the 1/N C term in 
the SU(N C ) Fierz identity, 

Y,(Xi\ a X 2 ){Y l \ a Y 2 ) = (X 1 Y 2 )(Y l X 2 ) - ^(X 1 X 2 )(Y 1 Y 2 ), (3.14) 

a c 

where AQ, Y{ are strings of generator matrices A ai . Such manifest 1/N C terms are not generated 
by strings of structure constants. (In fact, 1/N C terms are generally suppressed when / s are 
present, because l/N c factors come from explicitly projecting out the U(l) factor in U(N C ) = 
SU(N C ) x U(l), while structure constants accomplish this projection automatically.) Because of 
this fact, the freedom to include structure constants in a color decomposition does not seem to lead 
to any simpler representation of the subleading rif terms in eq. ( [3.12|) , and we shall be content to 
simplify the no-n/ termsQ 

The new color decomposition for one-loop amplitudes with an external qq pair is 



^ l0 ° P fe, 92, 93, ...,3n) = 9 n 



Y^ Yl ( AX2A<lCT3 • • • A^A* 1 )^ 1 ^"^ 1 • • • Fa ° n )x lX2 

p=2 <reS n - 2 

x^'W (lq, cjp+i, . . . , a n , 2 q , 0-3, . . . , a p ) (3.15) 

n-X 

+ ^E E Gv^(as,...,a n )A^(l q ,2 q ;as,...,a n ) 

j=i a-eS n -2/S n -j 

where for p = n the product of generators in the adjoint representation reduces to the identity, 
(F ■ ■ ■ F) XlX2 — ► S XlX2 ; the color structures Gr^ are defined in eq. ( 3.10 ); and 

A£{? ] (l q ,2 q ,Z,...,n) = A^ 2 \l q ,2 q ,^...,n), (3.16) 

A^ J 2] (l q ,2 q ;3,...,j + l;j + 2,j + 3,...,n) = (-1)^ ^ A^ l ' 2 \a{l qi 2„ 3, . . . ,n)) , 

aeCOP{a}{/3} 



2 Since the ty terms in the A n -j>i subamplitudes arise from projecting out a (7(1) factor, they can also be viewed as 
subamplitudes for a fictitious SU(Ni) x SU(N2) x U(l) theory containing two fermion representations, say (Ni, l)+i 
(the external qq pair) and (1, Na)+i (the fermion in the loop). This description makes it a bit more transparent how the 
color structure tr(A a3 ■ • ■ X aj+1 ) x (X aj + 2 ■ ■ ■ A"™)^ 1 multiplying these subamplitudes corresponds to their definition. 
But it does not lead to a distinct color decomposition. 
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with {a} and {/3} as in eq. ( |3.12| ). The new types of color structures appearing in eq. ( |3.15| ) are 
shown in fig. ||b. 

Most of the subamplitudes A n do not contribute at leading order in the large N c limit. (Only 



i-RJi 



those of the form A n (!«> • • • , 2„) do.) Therefore the technique of contracting with a suitable color 



structure at large N c cannot be used to prove eq. (|3.15[) . Instead we follow the first line of reasoning 
that we used at tree-level; that is, we substitute for the structure constants in terms of A matrices 
and collect terms. The result agrees with eq. (|3.9|), after substituting in eqs. ( |3.11 ) and ( ^.12[) for 
the leading and subleading subamplitudes and making use of the reflection identity ( 3.13[) . 

The contribution of one-loop amplitudes with an external quark-antiquark pair to NLO QCD 
cross sections parallels the n-gluon case, eq. ( |3.4j ). Carrying out the color-summed NLO interference 
terms for amplitudes color-decomposed as in eqs. ( |2.2| ) and ( 3.15J ), and omitting the nf terms, we 
obtain, 



^ A(qi,q2, 93,---, 9n)no-n f [A(qi,q 2 , g 3 ,---,9n 

2 y, Re ( At 



-Uf\ 



I NLO 



colors 



i tree 



(A 



1-loOpN 



colors 



2( 5 2 r- 1 Re^ 



n (n-2)! 



/no— nt 



(p) A tree. 



(3.17) 



2 i,j=l 



m-j'A 



*-qi £qi *3> 



R,[i], 



where 



- (p) 



colors 



Pi{\ a3 ■ ■ ■ X^jV 1 Pj{(X X2 X a3 



A 



A^A* 1 ).* 1 



qi 3p+li 



j 3ni Aqi .73) ■ ■ • j Jp) 



(pap + l . . . pa„ 



IXIX2J 



(3.18) 



with Pi the i permutation in S n -2, and Pj the j permutation in S n -2 a t fixed p. The partitions 

(v) 
p span a matrix- valued vector rh ), , which we give in appendix A for n = 4, 5, as required for cross- 



section computations up to 0{a*). For n = 5, we have checked that eq. ( |3.17| ), with mf- given in 
eq. ( [A. 16 ), agrees with the no-n/ terms of ref. [13]. The ny terms for n = 4, 5 have been computed 
in refs. |ll| and |D|, respectively. 

Finally, we note that by converting gluons into photons we can obtain multi-photon amplitudes at 
tree and loop level. In order to convert a gluon into a photon we make use of the embedding of U{\) in 
U(N C ). We replace the quark-gluon vertex factor gA" with the quark-photon vertex factor \^2eQ g I, 
where Q q is the quark electric charge and the factor \[2 is due to our choice of normalization of the 
A matrices. Since the identity matrix / commutes with the A matrices, all possible attachments of 
the photon to the quark line contribute to the same color structure. Thus subamplitudes involving a 
photon can be obtained by summing over permutations of gluon subamplitudes in which the photon 
assumes all possible places in the ordering. 

Photons do not couple to gluons. In a / -based decomposition this is apparent: When the 
index a corresponds to the identity matrix /, the U(N C ) generator F a in the adjoint representation 
vanishes. In this paper we presented new color decompositions for gluon tree amplitudes, eq. ( |2.1| ), 
and for the gluon-loop sector of one-loop amplitudes, eqs. ( ^ ) and ( |3.15| ); we still use the A-based 
decomposition for quark-gluon tree amplitudes, eq. (|2.2[), and for the quark-loop sector of one-loop 
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amplitudes. Thus it might seem that we can bring no new insight to the calculation of multi-photon 
amplitudes. However, that is not the case for one-loop amplitudes with an external qq pair: In the 
no-n/ terms in eq. ( 3.15| ) a photon can still attach to the external quark line, i.e. to the string of 



A matrices; however, it manifestly cannot attach to the string of F matrices. We immediately see 
that subamplitudes of type An (1$, 03, • • • , &n, 2<j), corresponding in eq. ( |3.15| ) to partitions with 
p = 2, do not contribute, because there is no A a matrix to replace by the identity. (The \ Xi matrices 



correspond to internal gluons and so should not be replaced.) In the A-based decomposition [2C], the 
vanishing contribution of these subamplitudes is realized only after summing over all permutations. 
Thus the new color decomposition we obtain for qqg ■ ■ ■ #7 one-loop amplitudes, by converting a 
gluon into a photon in eq. ( |3.15| ), and summing over all the possible attachments of the photon 
to the quark line, is more compact than the one obtained from the A-based decomposition. For 
n = 5 we verified that explicitly by computing the qqggj one- loop amplitude from eq. ( 3.15| ), and 



comparing it to (and checking its equivalence with) the same amplitude, but A-decomposed 

4 Discussion 



In ref. [|14| a new color decomposition was presented for the n-gluon tree amplitude; it was shown 
to be equivalent to the standard one up to n = 7 and was conjectured to be valid for arbitrary n. In 
the present paper we have proven this conjecture. We have also presented, and proven, analogous 
color decompositions for one-loop n-gluon amplitudes and one-loop amplitudes with an external 
quark- antiquark pair plus n — 2 gluons. 

The new color decompositions have the advantage of being written in terms of just the indepen- 
dent subamplitudes. This makes the decomposition of an amplitude more compact. It is particularly 
useful for computing the color-summed interference terms between tree amplitudes and one-loop am- 
plitudes, and the square of one-loop amplitudes, which are relevant respectively for NLO and NNLO 
calculations of jet production rates. The color-summed interference terms and squares are given in 
terms of color matrices; we have provided in the appendix explicit values required for cross-section 
computations up to 0{a A s ). 

We expect that structure-constant-based decompositions will also be useful for general multi- 
loop scattering amplitude calculations, once they have been derived. In a particular case they have 
already proven useful: Two-loop four-gluon amplitudes in N = 4 supersymmetric Yang-Mills theory 
were computed via their unitarity cuts and were presented in a trace-based formalism [^] , but they 
can be more compactly expressed in an f abc basis [22]. In the latter basis only two types of color 
structures appear; these are represented graphically in fig. ||a by the planar and nonplanar double 
box diagrams. The kinematic object multiplying each structure is simply the scalar integral (i.e., 
the 3 Feynman diagram) for the graph with the same topology. 

Some of this simplicity is undoubtedly due to N = 4 supersymmetry. However, group theory 
alone restricts the possible multi-loop color structures considerably. At one loop, we used the Jacobi 
identity to remove nontrivial trees from the set of n-gluon color structures. We can similarly remove 
all such trees from multi-loop n-gluon color structures. The remaining color factors are illustrated in 
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(a) 



(b) 






Figure 6: (a) The only color factors appearing in the two-loop 4-gluon amplitude in N = 4 super- Yang- 
Mills theory, (b) The types of color factors for multi-loop n-gluon amplitudes that remain after removing all 
nontrivial trees, (c) Example of an 'unwanted' 4-gluon planar color factor. 

fig. |6|b (for the two- loop case). In the one- loop case, this reduction exhausted the use of the Jacobi 
identity; in the multi-loop case this is no longer true. For example, in the two-loop four-gluon 
amplitude the nonplanar double box diagram in fig. ya could be removed in favor of the planar 
double box diagram in fig. |6|a and another planar diagram, shown in fig. |6|c. Yet this would be an 
unwise move, at least in the N = 4 supersymmetric case, in the sense that the nonplanar scalar 
double box integral would then have to be associated with planar color factors. Thus the Jacobi 
identity should be wielded more judiciously in the multi-loop case, and we leave this to future work. 

Finally, we remark that even if a general, multi-loop f a -based color decomposition is not yet 
available, the results in this paper can still be employed in the context of specific multi-loop gluonic 
calculations based on unitarity cuts [^j], [7L |22| . For example, the unitarity cuts of two-loop amplitudes 
are given in terms of tree-level and one-loop amplitudes. Color decomposing these amplitudes as 
in this paper leads directly, after sewing the color factors on each side of the cut, to an / -based 
decomposition for the cut of the two-loop amplitude. This decomposition is gauge invariant on 
the cut, by virtue of the gauge invariance of the tree and one-loop amplitudes. The only potential 
difficulty is that to reconstruct the amplitude from its cuts, one has to find a single function that 
consistently matches all the cuts, and the color factors may have to be manipulated at this point. 
Still, this seems likely to be a more efficient than a corresponding approach using a trace-based 
decomposition of tree and one-loop amplitudes, at least for computing terms that are subleading in 
the number of colors. 

At first sight it may seem like a step backward to abandon trace-based color decompositions of 
gluonic scattering amplitudes in favor of ones based on structure constants. However, we have seen 
in this paper that with an appropriate choice of color factors, one does not really have to give up 
any desirable properties of the trace-based formalism (such as gauge invariance) , at least at the tree 
and one-loop level. We look forward to multi-loop extensions of this work, and their application to 
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more precise jet rate computations. 
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A The color matrices 



We first give the color matrices appearing in the color-summed tree-level cross-section, defined in 



eq. ( J2.ll ). For n = 4 they are 



dij = C A {N C ) 



4 2 
2 4 



A\ 



tree 



(1,2,3,4) 

(1,3,2,4) 



(A.l) 



and for n = 5, 



c '.j 



C 5 (N C ) 



( 

4 
4 
2 
2 
\0 



8 4 



2 

4 
2 



o\ 

2 

2 
4 
4 
8/ 



i4 



tree 



/ (1,2,3,4,5) \ 
(1,2,4,3,5) 
(1,3,2,4,5) 
(1,3,4,2,5) 
(1,4,2,3,5) 

\ (1,4,3,2,5) ) 



(A.2) 



where C n {N c ) is defined in eq. ( |2.13|) . The above results agree with ref. jl7| , up to an overall factor 
due to the different choice of normalization of the matrices in the fundamental representation. The 
color matrix Cy for n > 5 shows the 0(1/N^) dependence outlined in eq. ( |2,9D . For n = 6, Cy 
is a 24 x 24 matrix. Rather than give the entire matrix, we note that its essential information is 
contained in its first row, i.e. in the matrix elements of (1,2,3,4,5,6) with the 24 permutations. 
Defining c (123 456)(i CT (2345)6) = C 6 (N C ) 7(cr(2345)), we have 



7(2345) = 
7(2354) 
7(2435) 
7(2453) 
7(2534) 
7(2543) 



16, 



4, 
4, 

0, 



7(3245) = 8, 7(4235) = 4, 

7(3254) = 4, 7(4253) = 2, 

7(3425) = 4, 7(4325) = 0, 

7(3452) = 2, 7(4352) = 0, 

7(3524) = 2, 7(4523) = 2 + a, 
7(3542) = 0, 7(4532) = a, 



7(5234) 
7(5243) 
7(5324) 
7(5342) 
7(5423) 
7(5432) 



2. 
0, 
0, 
a, 
a, 
a, 



(A.3) 



where a = 24/N^. The terms proportional to a agree with eq. (A. 10) of ref. p 

The color matrices for the tree- level one- loop interference, defined in eq. ( p.5|) , are for n 



C 5 (N C ) ■ 







dij = C A (N C 



-1 
-1 1 
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(1,2,3,4) 
(1,2,4,3) 
(1,3,2,4) 



(A.4) 
(A.5) 
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and for n 



C 6 (N C ) 



d ij =C 5 (N c 



[i] 



(4fi 





^2-2 





-2 





2 








a 


—a a a \ 




-2 


2 


-2 


2 





a 





a 


0a 







a 


-2 


—a 


2 


2 


-2 





a a 




a a 


a 








2 





-2 


2 


2 a 




a -2 








2 


—a 


a 





a 


2 2 




V a 


a 


a 


2 





a 


-2 





2 2/ 




I 1 - 1 





-1 





1 








6 


-b b b \ 




-1 


1 


-1 


1 





6 


6 


b 


6 







& 


-1 


-6 


1 


1 


-1 





6 6 




6 b 


6 








1 





-1 


1 


1 6 




b -1 








1 


-6 


6 





6 


1 1 




\ b 


& 


6 


1 





b 


-1 





1 1/ 


(1,2,3,4,5) 


\ 
















(1,2,3,5,4) 


















(1,2,4,3,5) 


















(1,2,4,5,3) 


















(1,2,5,3,4) 


















(1,2,5,4,3) 


















(1,3,2,4,5) 




5 














(1,3,2,5,4) 


















(1,3,4,2,5) 


















(1,3,5,2,4) 


















(1,4,2,3,5) 


















(1,4,3,2,5) 


/ 

















(A.6) 



(A.7) 



(A. 



V 

where a = 24/ N^, b = 2/N%. Although the matrices c™ and dL- have a relatively small dimension, 
they are not particularly sparse. We can improve the sparseness considerably by using the tree-level 
[7(1) decoupling identity, eq. ( |1.2j ), in every column where three identical entries appear, thus trading 
three terms for one. Of course, more than (n — 2)! different tree subamplitudes will now appear. 
Using also the tree-level reflection identity, we arrive at the formula 



5^A(1,... ,5) [A(l, 



■ 5)]« 



NLO 



(A.9) 



colors 



(n-l)!/2 

4 5 8 C 6 (iV c )Re £ 
i=l 



AF"lJlP + 2rAW\ 



N P . 



) 



,_L A tree /^.[l] , , ^L .[1/2] 

+ ^2^(24135) i{ L/A i +l Nc A i 



where (24135) • i is the permutation (24135) composed with the i permutation. In comparing this 
formula to eq. (11) of ref. P3| ], we see that the A^ terms can be removed, at the price of multiplying 
the remaining (l/iV 2 )-suppressed pure-glue (no nf) terms by a factor of six. This form for the NLO 



five-gluon cross-section has been noted independently by W. Kilgore [24]. 
For the NNLO terms in eq. ( |3.7| ) with n = 4 we have 

/ 2 + a 
Cij = C 6 (N C ) ■ a 

\ a 




(A.10) 



16 



d ij =C 5 (N c ) 



1 + 6 6 b 

b 1 + 6 b 
b 6 1 + 6 



(A.ll) 



/ l-b+d d d 

Cij =C 4 (N C ) • d l-b + d d 

\ d d l-b+d 

where a,b are defined below eq. ( |A.8| ). Equivalently, 
J2 lA^CM.M)! 2 



J_ _3_ 



colors 



9 S C 6 (N C ) 



£( 2 |^| 2 + 4^Re4 ] (4 1/21 )* +4^(1 - 6) |A^| 2 



i=l 

2 



A c 

3 



t=l c i=l i=l c t=l 



[1/2] 



» 



(A.12) 



(A.13) 



The color matrix-valued vector m^ for the tree-level one-loop interference for NLO processes 
with an external quark-antiquark pair, defined in eq. ( |3.18| ), is for n = 4, 



/ 



.(2,3,4) _ r (N . 



1 

1 o 
l/i i 



\ 



A c 2 V 1 1 
_L_ / e 1 + e 
V M l l + e e ; / 



(A.14) 



4 ree =( (1,2,3,4) (1,2,4,3) 



where e = 1/A 2 . For n = 5 we obtain, 



(A.15) 
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Ill 



(2,3,4,5) 
ij 



-C 5 (N C 



1 



N* 



f 

e + f 

e + f 
-l + 2e + / 
-l + 2e + / 

3e + / 



1 



b 

b 

b 



1 

/o 
1 

2 


1 
V i 



i 



1 
1 

2 




e 

e 

1 + e 

-1 + e 

e 

1 + e - 

e + f 

f 
-l + 2e + / 

3e + / 

e + f 

-l + 2e + f 



e 
e 
e 

1 + e 
1 + e 

-1 + e 



1 + e 
-1 + e 

e 

e 
1 + e 

e 

e + / 

1 + 2e + / 

/ 
e + / 

3e + / 
1 + 2e + / 





2 

1 

1 



1 

e 
1 + e 

e 

e 
-1 + e 
1 + e 



/ 

1 \ 

1 


2 
1 
0/ 



-1 + e 
1 + e 
1 + e 

e 

e 

e 



l + 2e + / 
3e + / 
e + / 

/ 
l + 2e + / 

e + f 



1 + e N 

e 

-1 + e 

1 + e 

e 

e y 

1 + 2e + / 

e + f 

3e + / 

1 + 2e + / 

/ 

e + f 



(A.16) 



3e + / \ 
-l + 2e + / 
-l + 2e + / 

e + f 

e + f 
f J 



A* 



(1, 2, 3, 4, 5) (1, 2, 3, 5, 4) (1, 2, 4, 3, 5) (1, 2, 4, 5, 3) (1, 2, 5, 3, 4) (1, 2, 5, 4, 3) , (A.17) 



where / = 1/N*. In both eqs. ( A.14J) and ( |A.16 ), for each value of p, the permutations Pj of the 
gluons in the loop amplitude are understood to appear in the same order as the tree amplitude 
permutations A^ ee . 
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